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We introduce a solvable model of a nonlinear double-barrier structure, described by a generalized 
effective-mass equation with a nonlinear coupling term. This model is interesting in its own right for 
possible new applications, as well as to help understand the combined effect of spatial correlations 
and nonlinearity on disordered systems. Although we specifically discuss the application of the 
model to electron transport in semiconductor devices, our results apply to other contexts, such as 
nonlinear optical phenomena. Our model consists of finite width barriers and nonlinearities are dealt 
with separately in the barriers and in the well, both with and without applied electric fields. We 
study a wide range of nonlinearity coupling values. When the nonlinear term is only present in the 
barriers, a sideband is observed in addition to the main resonance and, as a consequence, the current- 
voltage characteristics present two peaks at different voltages. In this case, the phenomenology 
remains basically the same through all the considered variation of the parameters. When the well 
is the component exhibiting nonlinearity, the results depend strongly on the nonlinear coefficient. 
For small values, the current-voltage characteristics exhibit lower and upper voltage cutoffs. This 
phenomenon is not present at higher nonlinearities, and eventually linear-like behavior is recovered. 
We complete this exhaustive study with an analysis of the effects of having simultaneously both kinds 
of nonlinear couplings. We conclude the paper with a summary of our results, their technological 
implications, and a prospective discussion of the consequences of our work for more complex systems. 



I. INTRODUCTION 

This decade is witnessing a rapid increase of the 
amount of effort devoted to the study of one-dimensional 
(ID) systems, for reasons both theoretical and technolog- 
ical. In particular, emphasis is being put on localization 
and derealization in non-spatially-periodic (quasiperi- 
odic or disordered) problems, which are relevant in very 
many contexts. From the fundamental viewpoint, re- 
search is being conducted in two different main direc- 
tions. On the one hand, the generality of localization 
phenomena in ID systems, has recently been questioned 
in a number of papers,tr0 which have shown that the 
introduction of short-range correlation in the distribu- 
tion of inhomogeneities leads to the formation of bands 
of extended states. This unexpected phenomenon can 
be understood in terms of the structure of the transmis- 
sion coefficient for different segments of the considered 
system.B On the other hand, nonlinear effects have also 
been proposed to counteract the localizing influence of 
disorder on waves of any nature (see, e.g., the reviews 
in Ref. ^ and references therein). Unfortunately, the re- 
sults of the interplay between disorder and nonlinearity 
depend to a certain extent on the specific model consid- 
ered and a general theory of this problem is still lacking 
(see, for instance, Refs. |j and || and references therein). 

Following along the above line of research, in this pa- 



per we concern ourselves with the study of a nonlinear 
double-barrier resonant structure. Recently, several au- 
thors have reported research on related devices: For in- 
stance, a mean field analysis of multiple resonant tunnel- 
ing exhibiting chaotic behavior has been carried out by 
Presilla, Jona-Lasinio, and Capasso;Q in a different con- 
text, directional couplers including nonlinear elements 
have been investigated as optical switches.13 The model 
we propose shares some of the characteristics of these: 
Specifically, it is described by a Schrodinger equation 
where an effective nonlinearity is included and, being 
formally close to effective nonlinear Schrodinger equa- 
tions arising in other problems, it is also very general. 
However, we are going to be mostly interested in its ap- 
plication to resonant tunneling (RT) through semicon- 
ductor heterostructures. This phenomenon, which takes 
place in linear double-barrier structures (DBS), make 
these systems very promising candidates for a new gener- 
ation of ultra-high speed electronic devices: For instance, 
a GaAs-Gai-^Al^As DBS operating at THz frequencies 
has already been reported in the literature.H 

The basic reason for RT to arise in DBS is a quan- 
tum phenomenon whose fundamental characteristics are 
by now well understood: There exists a dramatic in- 
crease of the electron transmitivity whenever the energy 
of the incident electron is close to onc,-pf the unoccu- 
pied quasi-bound-states inside the wells In practice, a 
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bias voltage is applied to shift the energy of this quasi- 
bound-state of nonzero width so that its center matches 
the Fermi level. Consequently, the j — V characteristics 
present negative differential resistance (NDR). In actual 
samples, however, the situation is much more complex 
than this simple picture. This is so even in good-quality 
heterostructures, when scattering by dislocations or sur- 
face roughness is negligible. In particular, inelastic scat- 
tering is always present in real devices. Examples of in- 
elastic scattering events are electron-lattice and electron- 
electron interactions, in which the energy of the tunneling 
electron changes and the phase memory is lost. Strong lo- 
cal coupling between electronic and vibrational degrees of 
freedom (i.e., when electrons propagate in a deformable 
charged medium) leads to a self-induced attractive force 
which we will account for by means of an effective non- 
linearity as discussed below. Moreover, Hartree-type re- 
pulsive interactions between electrons can be viewed as 
a self-induced repulsive force within an effective medium 
framework. There are currently available very many ma- 
terials with different optical or electrical nonlinear prop- 
erties, and then it may be possible to purposely build 
nonlinear DBS structures like those we work with. In 
view of this, the question arises as to whether such non- 
linear devices will have characteristics of interest for ap- 
plications. In any case, it seems that the (intentional or 
not) introduction of nonlincarity is going to have non- 
trivial consequences on the linear phenomenon of RT, 
and our purpose here is to study how RT is so modified 
as well as to analyze its measurable consequences. 

Aside from the above technological and experimen- 
tal motivations, the research we summarize in this pa- 
per has further goals. We have recently studied how 
short-range spatial comaLations affect electron localiza- 
tion, both theoreticallyp'El (even including possible three- 
dimensional effects, see Ref. [lCjl and from the viewpoint 
of real nanoelectronic devicesEJ In this paper, what we 
also intend to address is the question of how nonlincar- 
ity modifies our conclusions about linear correlated dis- 
ordered systems. We have already posed this problem 
in a preliminary workjij where we considered concen- 
trated nonlinearities using point-like barriers described 
by (5-function potentials. Nevertheless, such a situation 
is rather academic and, in addition, we have shown, for 
the case of linear correlated systems, that the results 
change dramatically when the width of the barriers is 
not neglected. til Therefore, similar variations may (and 
should be expected to) occur in the nonlinear problem 
when the barrier widths are taken into account, and this 
has to be done if any prediction is to be relevant for ac- 
tual experiments or devices. 

The paper is organized as follows. In Sec. II we present 
our model, obtained by including a nonlinear coupling in 
a generalized effective-mass equation. We particularize 
it for a DBS but we insist that the model is quite general 
and applicable in different physical contexts. We discuss 
in detail the physics underlying our choice for the cou- 
pling. We sketch the exact solution of the model, as well 



as the way to obtain the transmission coefficient as a 
function of the nonlinear couplings and the applied volt- 
age V. For completeness, we include an appendix con- 
taining a brief discussion of the range of applicability of 
the equation and its connection with physical interpre- 
tations. Afterwards, Sec. Ill contains the main results 
and discussions of our analysis concerning the applica- 
tion to semiconductor heterostructures, namely electron 
transmission and j — V characteristics. We also analyze 
how the results are changed by the simultaneous pres- 
ence of both types of nonlinearity. Finally, Sec. IV con- 
cludes the paper with a brief survey of the results and 
some prospects on the application of the ideas we have 
discussed. 



II. MODEL 

A. Physical grounds and definitions 

As we announced in the introduction, to describe our 
model we have chosen to apply it to a specific system: 
A semiconductor DBS under an applied electric field. In 
the following, we will be using parameters correspond- 
ing to GaAs-Gai-^Al^As as typical values for the linear 
DBS, whereas the nonlinear terms will be taken at most 
as a few percents of the built-in potential. Larger nonlin- 
earities are not considered because they may invalidate 
from the beginning the approximations involved in our 
theoretical calculation and, besides, they may make im- 
possible to find any material with the desired properties, 
at least in the context of electron transport. Our choices 
are thus as follows. The thickness of the whole structure 
is L and the thickness of the well is d. The barriers are as- 
sumed to be of the same thickness (symmetric case) but 
as will be evident below this is not a restriction of our 
approach. The structure is embedded in a highly doped 
material acting as contact, so that the electric field is ap- 
plied only in the DBS. We focus on electron states close 
to the bandgap and thus we can neglect nonparabolic- 
ity effects hereafter. Then the one-band effective-mass 
framework is completely justified to get accurate results. 
For the sake of simplicity, we will further assume that the 
electron effective-mass m* and the dielectric constant are 
the same in both materials. This hypothesis is related to 
the fact that we are not interested in high quantitative 
accuracy, although we note that the spatial dependence 
of these parameters can be taken into account if neces- 
sary. 

Within this approach, the electron wave function is 
written as a product of a band-edge orbital with a slowly 
varying envelope-function. Therefore the envelope- 
function ip{z) satisfies a generalized effective-mass equa- 
tion (we use units such that energies are measured in 
effective Rydberg, Ry*, and lengths in effective Bohr ra- 
dius, a*, being 1 Ry* = 5.5 meV and la* = 100 A in 
GaAs) given by 



2 



- i> zz (z) + [V(z) - eFz] %l>(z) = E $(z), (I) 

where V(z) is the potential term which we discuss below 
and F is the electric field applied along the growth direc- 
tion. We note, in connection with the generality of our 
model, that equations similar to Eq. (|]) are used to de- 
scribe light or other electromagnetic wave propagation in 
dielectric superlattices, in an approximation where.-pnly 
the scalar nature of the waves is taken into account .0 We 
now specify our model by choosing what is the potential 
term V(z). In order to do that, let us first consider the 
physics we are trying to represent with this term. The 
DBS can be regarded as an effective medium which re- 
acts to the presence of the tunneling electron, leading to 
a feedback mechanism by which inelastic scattering pro- 
cesses change the RT characteristics of the device. It thus 
follows that V{z) must contain nonlinear terms if it is to 
summarize the medium reaction which comes from the 
electron-electron and electron-lattice interactions. The 
simplest candidate to contain this feedback process is the 
charge density of the electron, which is proportional to 
^(z)! 2 - In our model, we neglect higher order contribu- 
tions and postulate that the potential in Eq. ([j]) has the 
form 

V(z) = V Q { [1 + &\^(z)\ 2 } X b{z) + mz)\ 2 Xw (z)} , (2) 

where Vb is the conduction band-offset at the interfaces, 
and Xb( z ) an d Xw{z) are respectively the characteristic 
functions of the barriers and the well, 

( 1, if 0< 2 < (L - d)/2, 
Xb {z)= I 1, if (L + d)/2 < z < L, (3a) 
[ 0, otherwise, 

. f 1, if (L - d)/2 <z<(L + d)/2, 
Xw(z) - | otherwise. ^ DJ 

and all the nonlinear physics is contained in the coeffi- 
cients a and (3 which we discuss below. 

There are two factors that configure the medium re- 
sponse to the tunneling electron. First, it goes without 
saying that there are repulsive electron-electron Coulomb 
interactions, which should enter the effective potential 
with a positive term proportional to the charge, i.e., the 
energy is increased by local charge accumulations, lead- 
ing to a positive sign for a and (3. On the other hand, 
in polar semiconductors, the electron polarizes the sur- 
rounding medium creating a local, positive charge den- 
sity. Hence the electron reacts to this polarization and 
experiences an attractive potential, which implies a and 
(3 negative. This happens, for instance, in the polaron 
problem in the weak coupling limit, which becomes valid 
in most semiconductors, and where it xan be seen that 
the lowest band energy state decreasesJlj It is then clear 
that in principle any sign would be equally possible for 
the coefficients if a and (3 are to represent the combined 



action of the polarization of the lattice along with repul- 
sive electron-electron interactions. Intuitively, however, 
it is most realistic to think that a will be negative, be- 
cause a positive nonlinear interaction would arise from 
negative charge accumulation in the barriers, which is not 
likely to occur. However, as far as the well is concerned, 
this is not so, because charge does tend to accumulate 
between the two barriers. This is the case considered, for 
instance, in the works of Presilla et aln where they in- 
troduce a term in the Schrodinger equation proportional 
to the total charge in the well (the integral of the square 
of the wavefunction) . Notice as we are going to apply 
a field to the DBS, high charge values will not build up 
between the two barriers, at least for low doping levels. 
Therefore, we will assume that lattice polarization effects 
are stronger than electron-electron interactions inside the 
well, thus leading to a negative (3. On the other hand, 
we discuss below mathematical reasons imposing that a 
and (3 have to be negative as expected, allowing for it 
to be positive only if they are small. This will be shown 
there to stem from the fact that we are studying a bound- 
ary value problem for an ordinary differential equation, 
and therefore the case of positive couplings is restricted 
to a fully dynamical study, which will be the subject of 
further work. til 



B. Analytical results 

We now work starting from Eq. (|l|) with the definition 
in Eq. (||) to cast our equations in a more tractable form. 
For simplicity, and because we are interested in intrinsic 
DBS features, we consider that the contacts in which the 
structure is embedded behave linearly. Therefore, the 
solution of Eq. (Q) is a linear combination of traveling 
waves. As usual in scattering problems, we assume an 
electron incident from the left and define the reflection, 
r, and transmission, t, amplitudes by the relationships 
( A ( e ik z + re -ik z\ z<0 
= 1 Ate*" z>L, (4) 

where fcp = E,k 2 L = E + eFL, and A is the incident wave 
amplitude. Now we define ip(z) = Atfi(z), a — <S|A| 2 , 
and (3 — f3\A\ 2 . Notice that a and (3 are dimensionless 
parameters and that they depend on the amplitude of 
the incoming wave, which will be relevant later. Using 
Eq. (|l|) we get 

-<j> zz {z) + [V{z)-eFz-E]<j>{z) = Q. (5) 

To solve the scattering problem in the DBS we develop 
a similar approach to that given in Ref. Ea. Since <p{z) 
is a complex function, we write <p(z) = q\z) exp[i^(z)], 
where q{z) and 7(2) are real functions. Inserting this 
factorization in Eq. (J5J) we have 7z(z) = q~ 2 (z) and 

- q zz (z) + —— + [VoXb(z) - eFz - E] q{z) + 
q A {z) 

+V [axb(z) + I3x w {z)]q 3 {z) = 0. (6) 
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This nonlinear differential equation must be supple- 
mented by appropriate boundary conditions. However, 
using Eq. this problem can be converted into a ini- 
tial conditions equation. In fact, it is straightforward to 
prove that 



q{L) = k L 1/ \ q z (L) = 0, 

and that the transmission coefficient is given by 

= 4fc og 2 (0) 

T 1 + 2M 2 (0) + fc o V(0) + 9 2 (0)<z 2 2 (0) ' 



(7) 



(8) 



Hence, we can integrate numerically (^|) with initial con- 
ditions (?]) backwards, from z = L up to z — 0, to obtain 
q(0) and q z (0), thus computing the transmission coeffi- 
cient for given nonlinear couplings a and /3, incoming 
energy E and applied voltage V — FL. 

Once the transmission coefficient has been computed, 
and recalling that contacts are linear media, the tunnel- 
ing current density at a given temperature T for the DBS 
can be calculated within the stationary-state model from 



m*ekBT 
2n 2 h 3 



f>QQ 

/ t(E, V)N(E, V) dE, 
Jo 



(9a) 



where N(E, V) accounts for the occupation of states to 
both sides of the device, according to the Fermi distribu- 
tion function, and it is given by 



N(E, V) = In 



1 + cxp[(£ F - E)/k B T] 
1 + exp[(E F -E-eV) /k B T] 

where ks is the Boltzmann constant. 



(9b) 
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III. RESULTS AND DISCUSSIONS 



In our calculations wc have considered a GaAs- 
Gao.65Alo.35As double-barrier structure with L = 3d = 
150 A. The conduction-band offset is Vq — 250 meV. In 
the absence of applied electric field and nonlinearities, 
there exist a single, very narrow resonance with r ~ 1 
below the top of the barrier, with an energy of 80.7 meV, 
and hence the well supports a single quasi-bound state. 
When voltage is applied, the energy of the quasi-bound 
state level is lowered and a strong enhancement of the 
current arises whenever the Fermi level matches this res- 
onance, thus leading to the well-known RT phenomenon. 
As we will see, it is enough to have a small amount of 
nonlinearity to have this picture changed dramatically. 
In this respect, before entering our report, we want to 
stress that the coefficients a and (3 that we have defined 
depend not only on the intrinsic characteristics of the ma- 
terial but also on the energy of the incoming electrons, 
as they both include the factor \A\ 2 . Therefore, it makes 
sense to study different values for the two coefficients, be- 
cause even for the same device transport properties can 
be very different for different incoming electrons. 




FIG. 1. Transmission coefficient r as a function of the elec- 
tron energy at zero bias for (a) a = —0.001, (b) —0.01, and 
(c) —0.1. For comparison, dashed line indicates the result for 
a = 0. 
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A. Nonlinear barriers 

We will begin by discussing the effects of having nonlin- 
earity in the barriers only (/3 = 0), and our starting point 
will be to report on zero field transmission properties. We 
will consider three values for a, a = —0.001, —0.01, —0.1, 
as representatives of very different orders of magnitude. 
Hereafter, for brevity of language we will call large values 
of the nonlinearity to those large in absolute value. For 
the purposes of the present paper we will regard values 
larger than those quoted above as not physically real- 
izable, although in principle materials with such exotic 
properties might be found; smaller values behave mostly 
like the linear system. Figure |l| shows the transmission 
coefficient as a function of the incoming energy for these 
three different values of the coefficient a, when the DBS 
is placed under zero bias. It is seen that for the two 
lower values there exist two peaks, but the smallest one 
is only well differentiated for a = —0.01. The main peak 
is centered at ~ 75meV (close to the linear case peak) 
and the sideband is centered at ~ 40meV for a = —0.01. 
Compared to the purely linear case, the main peak not 
only shifts to lower energies but also broadens, in a simi- 
lar fashion to what happens in RT when inelastic effects 
arise, and this shifting and broadening is larger when a 
increases. Upon further increasing a, the sideband moves 
below zero, whereas the main peak reduces appreciably, 
as shown in Fig. |l|(c). We thus see that the effect of 
strong nonlinearity becomes simply to move down and 
broaden the main peak, and therefore the existence of 
two peaks belongs only to a limited range of a. 

In order to gain insight into this phenomenon, we can 
rewrite Eq. (Q) as follows 

- 4>zz(z) + V eff (z, E) <f>(z) = E cf>{z), (10a) 

where we have defined an effective potential as follows 

V eff (z,E) = V oX b(z)[l + aq 2 {z)} - eFz. (10b) 

Thus ( |ic| ) is a Schrodinger-like equation for an effective 
potential due to nonlinearity plus the linear potential and 
the built-in potential of the DBS. This effective potential 
depends not only on z but also on the incoming energy 
E through the function q(z). Let us analyze its mean- 
ing in the zero field situation. Since the envelope func- 
tion changes under RT conditions, and also does q{z), it 
should be expected that V e ff(z, E) undergoes severe vari- 
ations whenever E is close to one of the RT peaks. This 
is indeed the (.'MSG , cLS shown in Fig. ||, where V e ff(z, E) is 
displayed at zero bias for a = —0.01 as a representative 
example of all the studied values. Notice that nonlinear 
effects have negligible effects on the shape of the effec- 
tive potential in the right barrier, besides a slight band 
bending at the interface z = (L + d)/2 at low energies. 
However, the potential in the left barrier region differs 
significantly from the original square-barrier shape. Out 
of energy resonances, the effective barrier height at z = 



is lower than Vq , whereas at resonances it takes the value 
~ Vq. This is clearly related to the ability of the electrons 
to go through the structure for the resonant energies and 
their subsequent low probability of remaining inside the 
barrier. Hence the effective potential presents two local 
maxima in the plane z = as a function of the incom- 
ing energy E, matching the values of the main resonance 
and the sideband above discussed. From Fig. ||it is clear 
that at the energy of the main resonance the effective 
potential is quite similar to the built-in DBS potential, 
just producing a small shift of the quasi-bound-state in 
comparison to the linear RT process. Concerning the 
sideband, the effective potential presents a deep mini- 
mum at the interface z = (L — d)/2, which originates 
another quasi-bound-state, thus explaining the origin of 
this lower RT peak. Indeed, the fact that this added 
well is responsible for the peak is confirmed by observa- 
tions of the effective potential when fields are applied, 
as discussed below. Another reason supporting this hy- 
pothesis is that for the smallest a the sideband is located 
at a lower energy, which makes sense since in that case 
the depth of the extra well is smaller. We thus provide 
a complete coherent picture of the zero field tunneling 
phenomenology. 

Effective potential (meV) ^^^^ \ JjB^&&^ 




FIG. 2. Effective potential V e // as a function of z and the 
incoming electron energy E at zero bias, for a = —0.01. 

We now consider the effect of a bias imposed on the 
DBS and show an example of the general behavior for 
all the considered values of a in Fig. 0. The transmission 
peaks are shifted to smaller values of energy, leading first, 
in the cases when there is a sideband, to the suppression 
of this subsidiary peak, as can be seen by comparing Figs. 
U (a) and (b). In this respect, we have to say that such 
shift to lower energies is similar to that found in linear 
(coherent) RT. Now we can conclude the interpretation 
of our findings in terms of the effective potential as an- 
nounced. By reasoning in the same way as above but 
taking the field into account, we have concluded that the 
fact that the resonance responsible for the peak shifts to- 
wards E = is due to the influence of the applied field, 
which increases the depth of the secondary well; the sec- 
ond resonance goes subsequently below zero and ceases 
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to be seen in the transmission coefficient. Afterwards, 
for higher fields, even the main peak is suppressed by the 
same mechanism, although if even higher fields are ap- 
plied, new levels are bound by such a deeper well. We do 
not go further into the details of those as we are inter- 
ested in the small fields suitable for our approximations 
and for application to actual devices. 
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FIG. 3. Transmission coefficient r as a function of the elec- 
tron energy for a = —0.01 with (a) V = 0.05 volts and (b) 
0.10 volts. 



We are now in a position where we are able to comment 
on the j — V characteristics, computed from (||) . We have 
set two different temperatures (77 K and room tempera- 
ture) and compared these curves with those obtained in 
linear RT. The Fermi energy was Ef — 27.5 meV. We 
discuss the lower a cases first, shown in Fig. ||. We have 
not plotted the curve for a = —0.001 as it turns out to 
be almost indistinguishable from the linear one, which 
always shows a single NDR peak. When a = —0.01, at 
T= 77 K nonlinearity causes the occurrence of a second 
peak at a lower voltage, clearly related to the sideband 



Applied voltage (Volts) 



FIG. 4. Computed j — V characteristics for Ef = 27.5 meV, 
a = —0.01 at (a) T — 77 K and (b) room temperature, and 
a = —0.1 at (c) T — 77 K and (d) room temperature. For 
comparison, dashed line indicates the result for a = 0. 
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On increasing temperature up to room values, both 
peaks merge into a single one because of the broadening 
of the Fermi-Dirac distribution function. It is important 
to notice that the j— V characteristics collect all the main 
features found in the transmission coefficient, namely a 
downwards shifting of the main RT peak, its correspond- 
ing broadening due to inelastic effects and the appearance 
of the sideband, as shown in Fig. ||(a) and (b) as a smaller 
NDR signature at V = 0.05 volts. As for the results for 
a = —0.1, they exhibit related features although some- 
what less prominent and shifted towards lower voltages 
[see Fig |](c) and (d)] . The disappearance of the sideband 
for this nonlinearity value leads to a structure of the j — V 
curve different than the previous one, which consists of a 
plateau starting at low voltage values. This is connected 
to the nonzero value of the transmission coefficient [see 
Fig. |l|(c)] for all energies, which allows for transmission 
even in non-RT conditions. The end of the plateau is an 
NDR interval which is associated to the broad peak orig- 
inated by the linear resonance. In addition, the current 
is twice as large as for the previous case. Following the 
plateau, the current drops down but it does not reach 
zero values, rather it keeps at a value about half that of 
the plateau. This is even more clear from the room tem- 
perature plot, which shows practically no features, i.e., 
loosely speaking, the current tends to be constant. The 
current values the DBS supports are clearly larger than 
in the linear and weakly nonlinear cases. 



B. Nonlinear well 



Let us now turn our attention to the case when a = 0, 
i.e., only the well in the DBS is nonlinear, which turns 
to be out very different from the previously discussed 
nonlinear barrier DBS. Results for the transmission coef- 
ficient at zero field are shown in Fig. ||. From those plots, 
it can be seen that in general the effect of (3 is more dra- 
matic than that of a, and that the case when f3 — —0.001 
has not much to do with the other two values. For the 
smallest nonlinearity, we see two peaks (as in the interme- 
diate a case), but the main one has departed largely from 
the linear resonance. This is related to the fact that (3 is 
negative, and therefore, accumulation of charge inside the 
well makes it deeper, thus moving the resonant quasilevel 
to lower energies. Increasing (3 will eventually locate the 
resonance below the zero energy and thus the appear- 
ance of the plots for (3 = —0.01 and (3 — —0.1. The less 
pronounced peaks appearing (meaning, for (3 = —0.001, 
the one playing the role of the sideband) have to do with 
the deformation of the bottom of the well by the spatial 
distribution of the charge. This can again be understood 
from the same effective potential ideas discussed in the 
previous subsection: By inspecting Fig. [| where the ef- 
fective potential is shown for the intermediate and large 
(3 values, it is clear that tunneling of electrons inside the 
well changes its structure dramatically, at first even split- 



ting it into two [Fig. ||(a)] or several [Fig. |6|(b)] wells. 
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FIG. 5. Transmission coefficient r as a function of the elec- 
tron energy for (a) (3 = -0.001, (b) -0.01, and (c) -0.1. The 
inset shows an enlarged view of the transmission coefficient 
close to 45meV. For comparison, dashed line indicates the 
result for /3 = 0. 
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The main point in that picture is that for the largest 
nonlinearity, the depth of the well is much larger than in 
the linear case, and therefore the main resonance is now 
located at negative energies, this being the reason for the 
disappearance of the main peak in the transmission coef- 
ficient. Besides, the complicated structure of the bottom 
of the well in [Fig. ^(b)] helps understand why the peak 
appearing in the transmission coefficient exhibits struc- 
ture; coupling between the subwells induces splitting of 
the quasi-level into several subresonances responsible for 
this finest features of the transmission coefficient [see the 
inset of Fig. ||(c)]. 



(a) 




value it is seen that one is entering from above. This 
will have consequences on the j — V curves which will be 
discussed below. 



(a) 
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(b) 
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FIG. 6. Effective potential V e ff as a function of z and the 
incoming electron energy E at zero bias for (a) f3 — —0.001 
and (b) -0.1. 



When bias is applied to the slightly nonlinear well 
(f3 = —0.001), the structure of the transmission coeffi- 
cient changes very smoothly, see Fig. 0, much as in the 
nonlinear barrier case. Due to the increased presence 
of electrons in the well induced by the field, the reso- 
nant level goes down and so does the transmission co- 
efficient peak. Remarkably, the intermediate nonlinear- 
ity (/3 = —0.01) is such that there is no peak at all in 
the transmission coefficient, although for the largest field 
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(e) 
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(f) 
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FIG. 7. Transmission coefficient r as a function of the elec- 
tron energy for j3 = —0.001 with (a) V = 0.05 volts and 
(b) 0.10 volts; (3 = -0.01 with (c) V = 0.05 volts and (d) 
0.10 volts; /3 = —0.1 with (e) V = 0.05 volts and (f) 0.10 volts. 

For f3 = —0.1, the field increases the inner structure 
of the peaks but for the rest they behave as in the other 
cases, moving downwards when increasing the field. In 
view of the effective potential structure, this probably has 
to do with the fact that when the field is applied the tilt- 
ing of the well induces the appearance of a positive energy 
quasi-level, whose energy decreases with increasing field 
and, subsequently, with increasing well depth. Larger 
fields destroy this structure, suppressing this quasi-level. 
The different behavior for the values of P considered re- 
flects in the j — V characteristics as well. We only show 
the results at T= 77 K because those at room tempera- 
ture are essentially the same, scaled down by a factor one 
half. In Fig. ||(a) it can be seen that the case /3 = —0.001 
is practically identical to the a = —0.01 case, which is 
reasonable in view that they both have the same trans- 
mission coefficient structure, albeit for different reasons. 
The main differences between those two devices are that 
for this nonlinear well one we are discussing, the first 
peak appears for much smaller values of the voltage, and 
that it shows a lot of structure coming from the inner 
wells induced in the quantum well. Upon increasing the 
value of /?, it is found that, interestingly, intermediate (3 
values suppress completely the current [Fig. ||(b)], a phe- 
nomenon that is evidently connected with the absence of 
resonances for that (3 value already reported. However, 
the case of large (3 shows a very unexpected feature: A 



series of noisy peaks appears when the nonlinear DBS 
starts conducing, which we tentatively associate to the 
structure of the well in this highly nonlinear case. For 
the rest, the DBS behaves qualitatively as in the linear 
case, as seen in Fig. ||(c). 



(a) 



Applied voltage (Volts) 




Applied voltage (Volts) 




Applied voltage (Volts) 



FIG. 8. Computed j— V characteristics for Ef = 27.5 meV, 
T = 77 K and (a) (3 = -0.001, (b) -0.01 , and (c) -0.1. Note 
that the scale in (b) is much smaller than in the other two 
ones. Results at room temperature at the same although 
scaled by one half. 
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C. Fully nonlinear DBS 

When both constituents of the DBS, i.e., the barriers 
and the well, are nonlinear, the number of possibilities 
is of course considerably augmented. It is not our aim 
to study thoroughly all possible combinations in this sec- 
tion, but rather, to present a brief idea of what can be 
said about this more complex devices. 



(a) 



25 50 75 100 

E (meV) 



that governs the behavior of the system. Thus, for 
(3 = —0.001, the structure of the transmission coeffi- 
cient shows two peaks until a — —0.1 which leads to a 
smoothing out of the peaks by almost suppression of the 
barriers. Indeed, the two peak structure is even clearer 
when (3 = a = —0.001 than in the other cases discussed 
so far [see Fig. ^(a)], and as consequence so is the two 
NDR peak feature of the j — V curve. The only differ- 
ence between a = —0.001 and a = —0.01 is that this last 
choice always gives rise to a somewhat more complicated 
curve [see an example in Fig. ||(b)]. On the other hand, 
when a is large, then it is this factor the one responsible 
for the behavior of the DBS, except when j3 is also large 
(thus compensating the decreasing of the barriers). In 
that shown in Fig. ||(c), it is possible to obtain 

a device capable of supporting a high current while still 
keeping a neatly marked NDR peak. We thus see that 
the desired characteristics of devices consisting only of 
one nonlinear component can be fine tuned by introduc- 
ing nonlinearity in the other component, or improved by 
suitably combining both. 
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FIG. 9. (a) Transmission coefficient r as a function of the 
electron energy for a = j3 = —0.001. (b) j — V curve for 
a = 10/3 = -0.01 and T = 77 K. (c) j -V curve for 
a = f3= -0.1 and T = 77 K. 

After an exploration of the main possibilities, we con- 
cluded that for small and intermediate a (barrier non- 
linearity), it is the nonlinearity of the well, (3, the one 



IV. CONCLUSIONS 

In this concluding section of the paper, we summa- 
rize briefly what are our main findings on the transport 
properties of nonlinear DBS. A first overall conclusion is 
that the behavior of such kind of devices changes quite a 
bit when the nonlinearity of the system changes. This is 
very important: Keeping in mind that the nonlinear co- 
efficients we have been handling include the contribution 
of the amplitude of the wavefunction, it is immediately 
understood that a device, built up from specific materi- 
als with specific and constant nonlinear characteristics, 
will change dramatically its response if the amplitude of 
the incoming wavefunction is also modified. This in turn 
can be related to the incoming current density. Clearly, 
a result like this paves the way to the search for devices 
with particular features arising from nonlinearity. We 
have to remind here that even if we have been discussing 
the outcome of our work in terms of electron transport, it 
can be straightforwardly translated to optical contexts, 
which gives further relevance to our results in view of the 
increasing importance of nonlinear optical phenomena for 
communications. In this respect, we want to stress that 
for that application the control of the nonlinearity of the 
devices by the incoming wave is even easier to realize. 

Being more specific, we have found that devices with 
nonlinear barriers and negative, i.e., self- attractive non- 
linearity, definitely show an increase of the current as 
compared to the linear case. This fact is not strange once 
it is realized that the role of nonlinearity in this case is to 
reduce the barriers. Roughly speaking, nonlinear barrier 
DBS begin transporting current for lower voltages; after- 
wards, they exhibit constant or slowly increasing char- 
acteristics for a range of voltages, which ends up in a 
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NDR region. The features of the j — V curves are of 
course smoothed when the temperature is rised up to 
room values, a conclusion which applies to all the non- 
linearities considered in this work. On the other hand, 
if the DBS is equipped with a nonlinear well, the device 
has more striking properties: In the lowest nonlinearity 
situation, it exhibits highly nonlinear properties, such 
as two NDR peaks; intermediate values of the nonlin- 
ear contribution will strongly suppress the current, and, 
finally, high nonlinearities lead to a quasi-linear current- 
voltage characteristic, with the RT peak being shifted to 
lower voltages and the current slightly increased. This is 
a peculiar property, and therefore we want to draw the 
readers' attention to the curious fact that nonlinear well 
DBS behave more differently from the linear case in the 
small nonlinearity situation. If we take also into account 
the suppression of the current for intermediate nonlin- 
ear values, it is evident that this kind of devices is the 
most appealing candidate for new applications. Finally, 
we have briefly pointed out some examples of combining 
both nonlinearities, and we have seen that it is possible 
to obtain clearer NDR peaks or higher conductivities (or 
both) by appropriately choosing the two components of 
the DBS to be nonlinear. 

In view of this report, we feel we can conjecture than 
complex structures consisting of a number of nonlinear 
DBS will also exhibit very peculiar properties. One can 
think, for instance, of a system built up of nonlinear well 
DBS. If the input to such device is of high amplitude, the 
first DBS components will behave quasi-linearly, letting 
the current pass through them. The deviation from unity 
of the transmission coefficient will make decrease the am- 
plitude, thus reaching the region for which the current is 
severely suppressed. This, in turn, would give rise to 
current in the remainder of the structure, because then 
it would behave more nonlinearly (as discussed above) 
allowing for transmission of the wavepacket. We would 
have thus designed a device which would transport cur- 
rent differently as a function of its length. Another possi- 
bly relevant remark relates to the mutual influence of dis- 
order and nonlinearity, and applies to the case of nonlin- 
ear barrier DBS. It is clear that if a superlattice made of 
linear DBS is built, the imperfections introduced during 
the growth will lead to a mismatch among the resonant 
levels of neighboring wells and as consequence to a loss of 
quantum coherence leading to poor transport properties. 
Nonlinear barrier DBS will be much more robust against 
growth imperfections, as they have a plateau-like current 
behavior for a range of voltages, which makes irrelevant 
the exact matching of quasi-levels. The conclusion is once 
more that nonlinearity would help transport even in the 
presence of disorder. In any event, a detailed study of 
more complex structures is necessary to clarify and put 
on firm grounds all these ideas. 

We finish by briefly pointing out another group of 
open questions which arise in view of our work. Fur- 
ther extensions of the present work to study nonlinear 
dynamical response of DBS on external ac bias would be 



of great interest to shed light on-related problems like 
bistability,E3 noise characteristics,!!!! and RT at far infra- 
red frequenciesta under the influence of inelastic scatter- 
ing channels as those described here. Besides that, the 
above mentioned highly nonlinear limit could also be in- 
teresting, for nonlinearity is always susceptible to give 
rise to new and unexpected features. If these new fea- 
tures were seen in our model it would be a very excit- 
ing development. If materials with suitable characteris- 
tics were found, and that is to be expected, experiments 
could be made to check the predictions: If the model were 
wrong, that would establish its range of validity, whereas 
if the predictions were correct, this work could pave the 
way to a new family of devices and applications. 



ACKNOWLEDGMENTS 

We are very thankful to Paco Padilla for useful discus- 
sions. A. S. acknowledges partial support from C.I.C. 
y T. (Spain) through project No. PB92-0248 and by 
the European Union Human Capital and Mobility Pro- 
gramme through contract ERBCHRXCT930413. 



APPENDIX: NONLINEARITY SIGN 

In this appendix, we give mathematical reasons why a 
and (3 (or a or 0) should be negative by using Eq. (||). 
We focus in the study of the problem with (3 = for sim- 
plicity, but similar considerations apply to the general 
case. With this assumption, the equation has the form 



q Z z{z) + 



1 



q*(z) 



+ f 1 (z)q(z) + f 2 (z)q 3 (z) = 1 (Al) 



where f\{z) and /^(z) are well-behaved functions. Let us 
now consider the discrete version of this equation, with 
the second derivative discretized in the usual way; denot- 
ing q n = q(z = nAz) and /,„ = fj(nAz), i — 1,2, with 



Az being the integration step, Eq. (Al) can be rewritten 
as 



(Az) 



1 



+ flnQn + hn<l n 



(A2) 



Notice that in this expression the sign of f2 n is the 
same as the sign of a and a. Let us now consider Eq. ( A2) 
for large q n ; this is general, because if q is small the 
term q~ 3 will make it grow quite quickly. In this limit, 



Eq. 



can be approximately replaced by 
Ag„_r = Aq n + (Azff 2n ql, 



(A3) 



where Aq n = q n — q n +i, and we have cast the equation in 
this fashion because it is to be integrated backwards. Re- 
calling the initial conditions (0) , if N is the total number 
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— 1/2 

of grid points, we have = Qn-i = k L > 0. There 
fore, we see that Agjv-i = 0, and 



Aq N _ 2 = {Az) 2 f 2n k 



-3/2 



(A4) 



We thus see that if f 2 n is positive (and hence a) the first 
increment is positive, and so are the subsequent ones, 
leading to a exponential divergence of q, whereas if f 2n 
and a are negative, the increment is negative, and q de- 
creases until the q~ 3 starts being relevant again. In this 
last case, it is possible that q reaches an equilibrium due 
to the balance of the two cubic terms, which was not for 
a > 0. This is in fact seen in the numerical integration 
of Eq. (^), where q rapidly diverges if a is positive, un- 
less, of course, a is positive but very small and/or the 
barrier (the region for a to influence q) is very narrow. 
In this last situation, however, the effect of a becomes 
negligible. 

Physically, this can be understood as follows. Seeing 
z as a time variable, Eq. (Al) can be regarded, loosely 
speaking, as an evolution equation for q, which is equiva- 
lent to the following interpretation: Electrons impinge on 
the barrier from outside, and begin to tunnel through it, 
their wavefunction being real and exponentially increas- 
ing or decreasing. If a is positive, then if there were any 
charge density in the BDS, this would become even more 
repulsive, and the wavefunction would diverge even faster 
(numerically one always see the exponentially growing 
part, of course). This instability is not present in the op- 
posite case, where a negative a helps the electron tunnel 
across the barrier. Note that this reasoning does not ap- 
ply to the full partial differential equation, which is the 
one we should deal with for that case by starting from 
the complete Schrodinger equation. 
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